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Abstract 

We study the persistence and propagation (or blocking) phenomena for a species 
in periodically hostile environments. The problem is described by a reaction-diffusion 
equation with zero Dirichlet boundary condition. We first derive the existence of a 
minimal nonnegative nontrivial stationary solution and study the large-time behavior 
of the solution of the initial boundary value problem. To the main goal, we then study 
a sequence of approximated problems in the whole space with reaction terms which are 
with very negative growth rates outside the domain under investigation. Finally, for 
a given unit vector, by using the information of the minimal speeds of approximated 
problems, we provide a simple geometric condition for the blocking of propagation 
and we derive the asymptotic behavior of the approximated pulsating travelling fronts. 
Moreover, for the case of constant diffusion matrix, we provide two conditions for which 
the limit of approximated minimal speeds is positive. 



1 Introduction and main results 



This paper is concerned with persistence and propagation phenomena for reaction-diffusion 
equations of the type 

u t -V -{A(x)Vu) = f{x,u) (f.l) 
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M-032-006-MY3. The second author is supported by the French "Agence Nationale de la Recherche" within 
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in W N or in some unbounded open subsets Q of R with zero Dirichlet boundary condition 
on dQ. Equations of the type (j 1.1ft arise especially in population dynamics and ecological 
models (see e.g. E2 EI]), where the nonnegative quantity u typically stands for the 
concentration of a species. 

Let us start with the case of the whole space M. N . The symmetric matrix field x i— > 
A(x) = i<ij<N is assumed to be of class C 1 ' a (lR iV ) with a > and uniformly positive 

definite: that is, there exists a positive constant > such that 

WxeR N , V£ =(&,...,£*) eR*, A£ ■ £ -.= ]T Mx)^ > p |£| 2 , (1.2) 

l<i,j<N 

where | • | denotes the Euclidean norm in Mr. We set R+ = [0, +oo). The nonlinear reaction 
term / : M N x R + — > M, (x,u) H- f(x,u) is assumed to be continuous, of class C 0,a with 
respect to x locally uniformly in a 6 R + , of class C 1 with respect to u, and §{(-,0) is of 
class C°' a CR N ). Furthermore, we assume that 



f(x, 0) = for all x G M N , 

there exists M > such that f(x, M) < for all x G 



(!-3) 



The functions Aij (for all 1 < i,j < N) and f(-,u) (for all u G R+) are assumed to be 
periodic in M N . Hereafter a function w is called periodic in M N if it satisfies 

+ k) = w(-) for all k G LiZ x ■ ■ • x L N Z, 

where Li, ■ ■ ■ ,Ln are some positive real numbers, which are fixed throughout this paper. 

If / fulfills the additional Fisher-KPP (for Kolmogorov, Petrovsky and Piskunov) [HI [20] 
assumption 

V x G M N , u i — y g(x, u) = — is decreasing with respect to u > 0, (1.4) 

u 

then the large-time behavior of the solutions of the Cauchy problem 

iH - V • (A(x)Vu) = f{x, u), t > 0, x G mN 
u(Q, x) = Uq(x), x G 



(1.5) 



is directly related to the sign of the principal periodic eigenvalue Ai of the linearized operator 
at (see [S]). This eigenvalue Ai is characterized by the existence of a (unique up to 
multiplication) periodic function (p G C 2,a (M N ), which satisfies 



-V • (A(x)Vtp) - C(x)<p = \ t <p inM N 

ip > in 



(1-6) 



where ((x) = ff(x,0) for all x G R N . The precise statement of what is known under the 
additional assumption (II. 4p will be recalled just after Proposition 11.11 below. 

Our first result, which is a preliminary step before the main purpose of the paper devoted 
to propagation phenomena in environments with hostile boundaries, is actually concerned 
with the existence of a minimal positive stationary solution p for problem (11.51) and with the 
large time behavior of the solutions u of (11. 5p . when / fulfills the assumption (11.31) alone. 
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Proposition 1.1 Assume that X% < and ( II. 3p . TTten there is a minimal periodic solu- 
tion p(x) of 

-V • (A(x)Vp) = f(x,p(x)) inR N , 

p > mR N , ^ ' ' 



in the sense that, for any solution q of (1.7), there holds q > p in R N . Furthermore, p < M 
in R N and, if Uq : R N — > [0, M] is uniformly continuous and not identically 0, then the 
solution u(t, x) of the Cauchy problem (\1.5§ is such that 



liminf u(t, x) > p(x) locally uniformly with respect to x G 

t— >+oo 



// one further assumes that u$ <p in R N , then u(t, x) — > p(x) as t — > +oo locally uniformly 
with respect to x G Mr. 

It is obvious to see that the solution p of (II. 7\i is not unique in general. Choose for 
instance A(x) = In (the identity matrix) and f(x,u) = sin(w) for all (x,u) G R N x R + : the 
function / satisfies (II. 3p with M = it, Ai = — 1 < 0, but any constant function p(x) = mir 
with m G N\{0} solves (II. 7p . On the other hand, if, in addition to ( II. 3p . the function / 
satisfies the assumption (ll.4p . then the solution p of (II. 7p is unique, see [5]. In particular, all 
solutions of ( II. 7p are necessarily periodic. Notice that, in the general case of assumption ( II. 3p 
alone, Proposition 11.11 still states the existence of a minimal periodic solution p of ( ll.Tjl in 
the class of all positive solutions q, which are not a priori assumed to be periodic. It is also 
known that, under hypotheses ( II. 3p and (II. 4p . the condition Ai < of the unstability of is 
a necessary condition for the existence of the solution p of ( II. 7p as well: if Ai > 0, then all 
bounded solutions u of ( II. 5p converge to as t — > +oo uniformly in R N , see [5]. On the other 
hand, under the assumptions ( II. 3p . ( II .4p and Ai < 0, for any non-zero bounded uniformly 
continuous u : R N — > R + , there holds u(t,x) — > p(x) as t — > +oo locally uniformly in 
x eR N (see [51 [02 ES])- We also refer to QUI [EE] for related results in the case of bounded 
domains with Dirichlet or Neumann boundary conditions, and to [U [7] for results with KPP 
nonlinearities and periodic or non-periodic coefficients in R N . Lastly, it is worth noticing 
that Proposition 11.11 and the aforementioned convergence results are different from what 
happens with other types of nonlinearities /, like combustion, bistable or even monostable 
nonlinearities which are degenerate at 0: in these cases, the large-time behavior of the 
solutions u of ( II. 5p strongly depends on some threshold parameters related to the size and/or 
the amplitude of u (see e.g. [II 121 [30l EH [38]). 

Remark 1.1 The assumption that u$ ranges in the interval [0, M] is made to guarantee the 
global existence and boundedness (from below by and from above by M) of the solutions u 
of the Cauchy problem (II. 5p . If / fulfills the KPP assumption (jl.4j) together with (ll.3p . or 
if f(x,s) < for all (x, s) G R N x [M, +oo), then it follows that the solution u exists for 
all t > and is globally bounded from below by and from above by max (M, ||wo||l°°(rjv)) , 
as long as uq is nonnegative and bounded. The same comment also holds for the Cauchy 
problem f 1 1 . 1 4 j) below with zero Dirichlet boundary condition on dQ. 
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As a matter of fact, in Proposition II. 1[ the negativity of Ai immediately implies that the 
positive periodic functions e ip are subsolutions of (jl.5p for e > small enough, where ip is a 
solution of ( II. 6p . It then follows from the above proposition and the results of Weinberger [32] 
that, for each unit vector e of M N , there is a positive real number c*(e) > (minimal speed) 
such that the following holds: for each c > c*(e), there is a pulsating travelling front 

u(t, x) = <p(x ■ e — ct, x) 

solving (11.51) and connecting to p, that is, the function (f> : IRxIR^ — > [0, M], (s, x) !->■ <p(s, x) 
is periodic in x, decreasing in s, and it satisfies <ft(— oo,x) = p(x) and <f)(+oo,x) = for all 
x G M N . Furthermore, such pulsating travelling fronts do not exist for any c < c*(e). We 
also refer to [2U EHl ESI EZ] for other results about pulsating travelling fronts in the whole 
space M , including other types of nonlinearities and the case of time-periodic media. 

Now, based on the previous results in Mr, we turn our attention to the main concern 
of this paper, namely the case when there are hostile periodic patches in the domain under 
consideration. We deal with persistence and propagation phenomena for reaction-diffusion 
equations of the type 

u t - V • (A(x)Vu) = F(x,u), i6fi, 
u(t, x) = 0, x G dfl, 

in an unbounded open set Q C M N which is assumed to be of class C 2,a (with a > 0) and 
periodic. The periodicity means that f2 — f2 + k for all k G L\$j x • • • x Ij-^'Ki. Furthermore, 
the fields A[x) and F(x,u) are assumed to be periodic with respect to x in Q, to have the 
same smoothness as before and to fulfill (11.2j) and ( II. 3p above, where x G M N is now replaced 
with j;6 0. In particular, assumption (ll.3p is now replaced with 

Fix, 0) = for all a; GO, 

y ' - 1.9 

there exists M > such that F(x, M) < for all x G tt. 

Throughout the paper, we denote 

C = On ([0,Li] x ••• x [0,L N ]) 

the cell of periodicity of Q. The zero Dirichlet boundary condition imposed on dQ means 
that the boundary is lethal for the species. Note that the unbounded periodic open set Q 
is not a priori assumed to be connected. The reason for that will become clear later, once 
the approximation procedure (jl.ISp below has been introduced. However, due to the global 
smoothness of dQ, the set Q has only a finite number of connected components relatively 
to the lattice L\L x ••• x LtvZ. That is, there is a finite number of connected compo- 
nents u\, . . . ,u m of Q such that u>i D (u>j + k) = for all 1 < % ^ j < m and for all 
k G L{Z x ■ • • x L N Z, and 

Q = |J Q h where ^ = (J + k. (1.10) 

l<i<m keLiZx---xL N Z 
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The sets are not uniquely denned, but the sets fit are unique (up to permutation), periodic, 
and Q { n Slj = for all 1 < z 7^ j < m. 

In the case of no-flux boundary conditions u(x) ■ (A(x)Vu(t,x)) = on dQ when Q is 
connected, much work have been devoted in the recent years to the study of propagation of 
pulsating fronts u(t, x) = <p(x ■ e — ct,x), where <f>(s, •) is periodic for all s G R and e is any 
unit vector, for various types of nonlinearities F, in straight infinite cylinders [H |3T] or in 
periodic domains (31 ESI E3 E2J E5]- In the case of KPP nonlinearities F, further properties 
of the minimal propagation speeds can be found in [HI E31 ESI ESI EB [23 E2J 135] . 

In this paper, we consider a larger class of reaction terms F, together with zero Dirichlet 
boundary condition. Let us first mention that, under the assumption that the equation ( 11. 8p 
is invariant in the direction x\ and under appropriate conditions on F, classical travelling 
fronts 

u(t, x) = (j)(xx - ct, x 2 , • • • , x n ) 

in straight infinite cylinders (in the £i-direction) with zero Dirichlet boundary condition 
are known to exist (see [24"| 154"] . including the case of some systems of equations). In this 
case, the profiles of these travelling fronts solve elliptic equations or systems. For pro- 
blem ( II. 8p in periodic domains, the reduction to elliptic equations does not hold anymore 
since the equation is not assumed to be invariant in any direction. Recently, existence 
results for problems of the type (jl.8p in connected two-dimensional periodically oscillating 
infinite cylinders with homogeneous isotropic diffusion (A(x) = I2) and KPP nonlinearities 
satisfying (jl.4p have been established, see [23]. In the present paper, the set Q is periodic 
in all variables xi, . . . , Xjy and the direction of propagation may be any unit vector e of R^. 
Actually, one of the novelties of this paper with respect to the previous literature is that the 
nature of propagation vs. blocking strongly depends on the direction e and on geometrical 
properties of the set Q itself. 

Let Xi^d denote the principal periodic eigenvalue of the linearized equation at in Q with 
zero Dirichlet boundary condition. That is, there exists a function ip e C 2 ' a (Q), which is 
periodic in Q and satisfies 

' -V • {A{x)Vip) - ({x)<p = \ liD ip 

<p = 

' tp > 

max ip > 0, 
n 

where ((x) = ^(x, 0) for all iGfl. If Q is connected, then ip > in f2 and ip is unique up to 
multiplication. Otherwise, in the general case, the function ip is unique up to multiplication 
in each set fli on which it is positive. More precisely, (p can be chosen to be positive on the 
(largest possible) set Q = [j ie j min ^i, where J TO j„ denotes the set of indices i G {1, . . . , m} 
for which the principal periodic eigenvalue Xi,q u d of the operator —V • (A(x)V) — ((x) in Qi 
with zero Dirichlet boundary condition on dQ{ is equal to Xi,d- That is, 

\ 1>D = min \\ n-D = Xi,n,.D for all i G I min . 



in fi, 
on dQ, 
in Q, 



[1.11) 
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The following theorem, which is analogue to Proposition II. 1[ is concerned with the ex- 
istence of a minimal nonnegative and non-trivial stationary solution of f ll.Sp in Q and the 
large-time behavior of the solutions of the associated initial boundary value problem, under 
the assumption that the steady state of (I1.8P is linearly strictly unstable. To do so, we 
introduce the set 

/_ = {ie{l,...,m}, A 1 ,n t n<oV (1.12) 



Theorem 1.2 Assume that X^u < 0, that is I ^ 0. Then there exists a minimal stationary 
periodic solution p(x) of 

-V • (A(x)Vp) = F(x,p(x)) inU, 

P = in dQ U \J i0 fli, (1.13) 

P > in Ui 6 /_A> 

in the sense that any bounded solution q of satisfies q > p in Q. Moreover, for any 

uniformly continuous function uq : Q — > [0, M] which is not identically 0, the solution u(t, x) 
of the initial boundary value problem 

u t - V ■ (A(x)Vu) = F(x,u), t > 0, x G H, 

u(t,x) = 0, t>0, xedtt, (1.14) 

u(0, x) = u (x), xeQ 



is such that 



liminf u(t, x) > p(x) (1-15) 

t— >+oc 



locally uniformly with respect to the points igO whose connected components intersect the 
support of uq. If one further assumes that uq < p in Q, then u(t,x) — > p(x) as t +oo in 
the same sense as above. 

As already emphasized, the periodic open set fl is not assumed to be connected, this is 
why the lower bound (11.151) or the convergence of u(t, x) to p{x) at large time can only hold in 
the (open) connected components C of the intersection of Q with the support of u (outside 
these components, the solution u(t, x) stays for all times t > 0). If such a connected 
component C is included in a set Qi with i G L, then Theorem 11.21 implies that u(t,x) is 
separated away from at large time, locally uniformly in C. However, fll . 15j) does not say 
anything about the behavior of u(t,x) when x G Ui^/_ ^ ip( x ) = there). Actually, for 
each Qi with i one has Ai^^d > and if F satisfies the additional assumption ( II. 4p 
in fli, then u(t, x) — > as t — )■ +oo uniformly in x G as follows from the same ideas as 
in [5]. 

The remaining part of this paper is concerned with the existence of pulsating fronts and 
the possibility of blocking phenomena for problem ( II. 8p with zero Dirichlet boundary condi- 
tion. The strategy, which is one of the main interests of the paper, consists in approximating 
the Dirichlet condition on dVL (and even in M Ar \fi) by reaction terms with very negative 
growth rates in R using the previous results and then passing to the singular limit in 
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the stationary solutions and in the pulsating travelling fronts as the growth rates converge 
to — oo in R-^O. This means that the quantity u lives in the whole space M N , but the space 
contains very bad regions. We will see that the location of the good vs. bad regions plays a 
crucial role in the dynamical behavior of the solutions. 

For this, let (/ n )neN be a sequence of real-valued functions defined in ~R N x R + such 
that each function f n : (x,u) (-> f n (x,u) is continuous, periodic with respect to x G M. N , 
of class C 0,a with respect to x G ~R N locally uniformly in u G R+, of class C 1 with respect 
to u with ff(-,0) G C°> a (R N ), and it satisfies (JT3J. Here we define N to be the set of all 
nonnegative integers. Furthermore, we assume that 

f n (x, u) = F(x, u) for all (x, u) G Q x R + and n G N, 

(f n ( x , u ))n£N is nonincreasing for all (x,u) G Q x R + , (1-16) 



g n (x,u) — oo as n — > +00 locally uniformly in (x,u) G (R x 



where 

g n (x, w) 



it u > 0, 

^ n -(x, 0) =: Cn(^) if u = 0. 



9m 

The last condition means that the death rate in the region R-^fi is very high, namely this 
region becomes more and more unfavorable for the species as n becomes larger and larger. 
Typical examples of such functions f n satisfying (11.91) and fll . 16j) are 

fn(x,U) = p n (x)u + /(it), 

where the function / : R + — > R is of class C 1 and satisfies /(0) = 0, f(M) < 0, and the 
functions p n : R^ — » R are periodic, nonpositive, of class C ,0,Q: (R Ar ), nonincreasing with 
respect to n, independent of n in Q, and p n — > —00 as n — > +00 locally uniformly in R iV \f2. 

For every n G N, let Ai in denote the principal periodic eigenvalue of the linearized operator 
at in R^. That is, there exists a (unique up to multiplication) periodic function ip n of 
class C 2 ' a (R N ), which satisfies 



-V • (A(x)Vip n ) - Cn(x)ifn = ^nPn in 



ip n > in R N . ( ^ 1 ' 17 ' ) 

We first establish the relationship between the principal eigenvalues Ai )n of (I1.17P and the 
principal eigenvalue Ai ^ of (II. lip , as well as the convergence of the minimal solutions p n 
of (jl.7p with nonlinearities f n to the minimal solution p of f ll .13)) . when Ai^ < 0. 

Theorem 1.3 Under the above notation, the sequence (Ai jn ) n£ N is nondecreasing and there 
holds Ai jTl Ai,£) as n — > +00. Furthermore, if X^d < 0, then the sequence (p n )neN of 
minimal solutions of 7\j with nonlinearities f n is nonincreasing and 



p n (x) — > Poo(x) as n — >■ +00 for all x G 
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where, up to a negligible set, p^o is nonnegative, periodic in R , p^o = in R the 
restriction of poo on Q is of class C 2 ' a (Q) and solves 

-V • (A(x)Vpoo) = F(x,poo) in ft, 
Poo = on dQ. 

Lastly, poo > p in Q, where p is given in Theorem M.B. 

We point out that, in general, the function poo is not identically equal to the solution p 
of (I1.13P in ft. However, it is well equal to p in ft if F fulfills (jl.4j) in ft. We refer to 
Remark 13.11 for more details. 

The last result is concerned with the asymptotic behavior as n — > +oo of the pulsating 
travelling fronts of the type <f) n (x ■ e — ct, x) connecting to p n (for problem (11 .ip in R N with 
nonlinearities f n ) and of their minimal speeds c*(e) > in any direction e (when Ai n < 0). 
The limit shall depend strongly on the direction e and blocking phenomena may occur in 
general. 

Theorem 1.4 Assume that A^d < and let e be any given unit vector ofR N . 

a) The sequence (c*(e)) n£ N is nonincreasing with limit c*(e) > 0. // all connected com- 
ponents C of Q are bounded in the direction e in the sense that 

sup | a; • e\ < +oo, (1-19) 

then c*(e) = 0. 

b) For any c > c*(e) with c > and for any sequence (c n ) nS N such that c n — )■ c asn — )■ +oo 
and c n > c*(e), the pulsating travelling fronts u n (t,x) = <p n (x ■ e — c n t,x) for (li.il) in W N 
with nonlinearity f n satisfy 



u n (t,x) -> 



u(t,x) in C\ and C% locally in R x O, 
m Lj oc (R x (R N \Q)) 



up to extraction of a subsequence, where u(t,x) = (p(x-e — ct,x) is a classical solution of (1 1 . £\) 
with Ut > inRxfi and <p(s, •) is periodic in Q for all s G R. Moreover, for any given i G 
one can shift in time the functions u n so that u(—oo, •) = and u(+oo, ■) > in fli. 

c) Assume here that A is constant. If there exist a unit vector e' ^ ±e and two real 
numbers a < b such that 

Q D {x G R N , a < x ■ e' < b}, (1.20) 

then c*(e) > 0. If there exist a unit vector e' , a point Xq G R and a real number r > such 
that e' is an eigenvector of A with e' ■ e ^ 0, and 

Q D {x G R N , d{x,x +Re') < r}, (1.21) 

where d denotes the Euclidean distance, then c*(e) > 0. 
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Theorem 11.41 provides a simple geometrical condition for the blocking of propagation, 
in a given direction e, in the presence of hostile periodic patches (by blocking, we mean 
that c*(e) y as tx y +00). Consequently, some quantitative estimates of the spreading 
speeds of the solutions u of the Cauchy problems (jl.5p with nonlinearities /„ can be de- 
rived. Indeed, for any compactly supported function u ^ 0, the solution u of (jl.5p with 
nonlinearity /„ spreads in the direction e with the spreading speed 



mm 



?e §iv-i i? . e> o £ • e 

in the sense that liminf t ^ +00 u(t, cte + x) > p n (x) locally uniformly in x if < c < to*(e), 
whereas \im t ^ +oa u(t,cte + x) = locally uniformly in x if c > u>*(e) (see [U H51 m 
particular, < w^(e) < c*(e). Hence, under condition (I1.19p . c*(e) = c*(— e) = and 
the solution u of (11. 5p with nonlinearity f n spreads as slowly as wanted in the directions ±e 
when n is large enough. In this case, since all connected components of Q are bounded in the 
direction e, pulsating fronts in the directions ±e for problem (11.81) in Q make no sense even 
if, under the notation of part b), the solutions u n can be shifted to converge to a non-trivial 
solution u of (11.81) in M x Q: what happens is that, in each connected component of Qj, u is 
just a time connection between and a non-trivial steady state. 

On the other hand, Theorem 11.41 also gives some simple geometrical conditions, of the 
types (ll.20p or (11.211) . for non-blocking in the directions ±e. These conditions mean that Vt 
contains a slab which is not orthogonal to e, or contains a cylinder in a direction which is 
not orthogonal to e. We do not know however if these conditions are optimal, even when A 
is constant. Lastly, Theorem 11.41 shows the existence of pulsating fronts for problem (11.81) 
in Q. Assume for instance that Q is connected, that is m = 1 under notation (11.101) . Then, 
there are pulsating traveling fronts, in the usual sense, in the direction e, connecting to a 
non-trivial periodic stationary solution of (jl.8p . Furthermore, if F is of the KPP type (11.41) 
in Q, the limiting state is unique and is then equal to the function p = given in Theo- 
rems [L2] and [L3] (see Remark 13. II below and the end of the proof of Theorem II .4p . However, 
Theorem 11.41 holds for general monostable functions F which may not be of the KPP type 
and it gives the first result about the existence of pulsating fronts with zero Dirichlet bound- 
ary condition in periodic domains (which may not be cylinders). 

Outline of the paper. Section [2] is devoted to the proof of Proposition 11.11 and Theo- 
rem 11.21 about the existence of minimal non-trivial stationary solutions p of problems (I1.7P 
and (I1.13P respectively, and about the large-time behavior of the solutions u of the Cauchy 
problems (ll.5p and (I1.14p . Section [3] is concerned with the proof of Theorem 11.31 and the 
relationship between the minimal solutions p n of problems (ll.7p with nonlinearities f n and 
the minimal solution p of problem (I1.13p . Lastly, in Section HI we do the proof of Theo- 
rem 11.41 and make clear the role of the geometrical condition (I1.19p in the blocking process 
as n — >■ +00. 
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2 Minimal stationary solutions and large-time beha- 
vior for the Cauchy problems (11.51) and (11.141) 



In the first part of this section, we first deal with the elliptic and parabolic problems (11. 7\i 
and (11.51) set in the whole space R N with the assumption (jl.3p on the nonlinearity /. Namely, 
we do the proof of Proposition 11.11 It is based on the elliptic and parabolic maximum prin- 
ciples and on the construction of suitable subsolutions. Since some parts of the proof are 
quite similar to some arguments used in [5] and [7], they will only be sketched. In the second 
part of this section, we will be concerned with the stationary and Cauchy problems (11. 1 3 [) 
and ( I1.14p posed in the set Q with zero Dirichlet boundary condition. That is, we will do the 
proof of Theorem II. 2\ which will itself be inspired by that of Proposition II. 1[ but additional 
difficulties arise. 

Proof of Proposition 11.11 Let <p be the unique periodic solution of (jl.6p such 
that max K iv(y9 = 1. Since the principal periodic eigenvalue Ai of ( II. 6p is assumed to be nega- 
tive, one can fix e G (0, M) so that f(x, s) > ((x) s + (Ai/2)s for all (x, s) e R N x [0,e Q ]. 
Now, for any e G (0, Eq], there holds 

- V ■ {A{x)V(e<p)) - f(x, eip) < -eV ■ (A(x)V(p) - ((x)e<p - = y^ < (2.1) 

for all x G R . In other words, the functions Sip are strict subsolutions of (jl.7p for all 
e G (0,e ]. 

Let now U be the solution of the Cauchy problem (jl.5p with initial datum Uq = ^oV 9 - 
Since < Uo < M and /(•, M) < in W N and since Uo is a subsolution of (II. 7ft . it follows 
that 

e ip(x) < U(t,x) < M for all (t,x) E R+ x R N 

and that U is nondecreasing with respect to t. Furthermore, by uniqueness for the Cauchy 
problem ( 11. 5p . U(t,-) is periodic in R N for each t > 0. From standard parabolic estimates, 
it follows then that 

U(t,x) — > p(x) as t — > +oo uniformly with respect to x G R , 

where p is a C 2,a (R N ) periodic solution of (I1.7P such that < e ip = U < p < M. 

Let us then show that p is the minimal positive solution of (I1.7P (in the class of all 
positive solutions of ( 11. 7ft . which are not a priori assumed to be periodic). Let q be any 
positive solution of ( 11. 7ft . Let \i,b{ v ,r),d denote the principal eigenvalue of the operator 

-V • (A(x)V) - C(z) 

in the open Euclidean ball B(y, R) of center y G R N and radius R > 0, with zero Dirichlet 
boundary condition on dB(y,R). For each point y G M. N and R > 0, the principal eigen- 
value \i t B(y,R),D is characterized by the existence of a function <p y ,R of class C 2,a (B(y, R)), 
solving 

-V • (A(x)Vcpy t R) - C(x)(fy tR = \i,B(y,R),D ¥y,R in B{y, R), 

Py,R > mB(y,R), 
¥ y ,R = on dB(y,R). 
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Up to normalization, one can assume that max B ^ R ^ (p VyR = 1, and the functions f ViR are 
then unique. As done in [7j, there holds 

^i,B( y ,R),D ->• Ai as R — > +00, 

uniformly with respect to y G R N . Since Ai < 0, one can then fix R > large enough 
so that \i,B( y ,R),D < Ai/2 for all y G R N . Thus, for each y G R N and e G (0, £q], the 
function £<p yt R satisfies 

-V ■ (A(x)V{£(p y>R )) - f(x,expy !R ) < -eV ■ (A(x)V(p v , R ) - C{x)e<p y , R - ^e(p ViR 

(\ Xl \ ( 2 - 2 ) 

= yM,B(y,R),D - y 

< 



in B(y,R). In other words, the functions ey^.R are strict subsolutions of (I1.7P in the 
balls B(y,R) for all e G (0, e ]. Now, fix y G R N and observe that min g ^ R ^ q > by 
continuity of q. It follows then that 

e* y : = sup je: G (0,e ], e(p V)R < q in 

is positive. We shall prove that e* = Eq. Assume not. Then < e* < Eq and E*ip VtR < q 
in B(y,R) with equality somewhere in B(y,R). Since q > and <p y>R = on dB(y,R), 
the functions E*<p y>R and g are equal somewhere at an interior point, in B(y, R). But E*<p ytR 
is a subsolution of ( 11. 7p . from (12. 2p . Since / is (at least) Lipschitz-continuous locally with 
respect to the second variable, uniformly in x, it follows from the strong elliptic maximum 
principle that E* y ip V)R = q in B(y,R), which is impossible since the inequality ( 12. 2 j) is strict. 
Therefore, e* = £q f° r & h 2/ G an d, in particular, 



But, by uniqueness of the principal eigenfunctions (p y>R and by periodicity of A and (, 
the function y !->■ (f y> R(y) is continuous and periodic in K.^. Since it is positive, one gets 
that minygiRjv <p y>R (y) > 0. Therefore, inf K iv g > 0. 
Define now 

e* = sup |e G (0, £r ], < <? m K^j, 

where we recall that ip is the unique periodic solution of (11.61) such that max K iv<£> = 1. 
Since q is bounded from below in the whole space M. N by a positive constant and since (p is 
bounded, one has e* > 0. Assume that e* < Eq. Then £*(p < q in and there exists a 
sequence (xA;)fceN in ffi^ such that 

£*(p(xk) — q(xk) — > as — > +00. 

By writing = x' fe + a/ fe ' with x' fc G L]Z x ■ ■ ■ L^Z and G [0, Li] x ••■ x [0, Ljv], it 
follows that the functions qk{x) = q(x + x' k ) converge, up to extraction of a subsequence, to 
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a solution q^ of (11. 7p such that e*(p < q^ in R N with equality somewhere in R N . As above, 
one concludes that e*ip = q^ in R N , which is impossible since e*(p is a strict subsolution 
of (II. 7p . from (12. ip . Therefore, e* = Eq, whence EQip < q in ~R N . The parabolic maximum 
principle implies that 

U(t,x) < q(x) for all (t,x) Gl+x R N , 

where we recall that U denotes the solution of (II. 5p with initial datum s <p. By passing to 
the limit as t — > +00, one gets that 



p(x) < q(x) for all x G 



Finally, let uq : R N — > [0, M] be a uniformly continuous function which is not identically 
equal to 0, and let u denote the solution of (11.51) with initial datum uq. The maximum 
principle implies that < u(t, x) < M for all (t, x) G R + x R N , and u(t, x) > for alH > 
and x G M^. With the same notation as above, there exists then e G (0, £0] such that 



e<fo,R <«(1,-) in 5(0, i?), 

where we recall that i? > was chosen so that \i : B( y ,R),D < Ai/2 for all y G M^. Let f be 
the solution of (II. 5p with initial datum 



v [x) 



E(p 0jR (x) if x G -8(0, R), 

if x G R N \B(0,R). 



Since < v < u{l, •) < M in E^, there holds 

< w(t, x) < w(t + 1, x) < M for all (t, x) G 1R+ x 



Furthermore, since f is a subsolution of ( 11. 7p because of (I2.2p and v = in M 7V \i?(0, i?), it 
follows from the maximum principle that v is nondecreasing with respect to t. Hence, from 
standard parabolic estimates, one gets that 



v (t, x) — » v^x) as t — > +00 locally uniformly in x G 



DiV 



where Voo is a solution of ( I1.7P satisfying t> < f 00 < M in M^. Notice in particular that v oo is 
positive in from the strong maximum principle, since v is nonnegative and not identically 
equal to 0. But the previous paragraphs yield then > p. Therefore, 

lim inf u(t, x) > p(x) 

t— !-+oo 

locally uniformly in x G R . Lastly, if Mo < P, then u(t,x) < p(x) for all (t, x) G R+ x M. , 
whence w(t, x) — >■ p(x) as t — > +00 locally uniformly in x G Mr. The proof of Proposition ll.il 
is thereby complete. □ 

Let us now turn to the proof of Theorem 11.21 Some of the ideas of the proof of Proposi- 
tion [1J] can be adapted. However, the case of problems ( I1.13P and ( I1.14p in Q is substantially 
more involved than the case of the whole space R N , mainly due to the fact that zero Dirichlet 
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boundary condition is imposed on dQ and the connected components of Q may be bounded 
or unbounded. 



Proof of Theorem 11.21 Remember that the sets fli given in f ll.lOp are all periodic and 
pairwise disjoint. We first work in each set for which \i,n h D < 0, that is i G We 
claim that, for each such index i 6 7_, there exists a periodic solution p^ G C 2,a (Qi) of the 
stationary problem 

-V-{A{x)Vpi) = F(x,pi(x)) in ^ 

Pi = ondtti, (2.3) 

Pi > in fli. 

Indeed, let cfi be the principal periodic eigenfunction of the operator —V • (A(x)V) — C(x) 
in Qi with zero Dirichlet boundary condition on dQ{. That is, the function tpi is periodic, of 
class C 2 ' a (Qi), and it solves 

-V • (A(x)Vtpi) - = Xi^D&i in fi», 

& = on dn t , (2.4) 

> in Qi. 

Up to normalization, one can assume that max^y^ = 1. Now, as in the proof of Propo- 
sition CCU since Xi^d < 0, there exists eo G (0, M] such that, for any e G (0, £q], the 
function e^i is a strict subsolution of ( \2.3L namely 

- V ■ (A(x)V(e&)) - F(x, e&(x)) < in Q { , (2.5) 

together with e^pi = on dQi and eifi > in But since the constant M is a supersolution 
of this problem, the solution m of the Cauchy problem 

(ui) t - V ■ (A(x)Vui) = F(x, m), t > 0, x G fi~, 

Ui {t,x) = 0, t > 0, x G dtli, (2.6) 

Ui(0,x) = e <pi(x), x G fi i; 

is such that £o^i(^) < x) < M for all (i, x) G (0, +oo) x and ttj is nondecreasing in i 
and periodic in x in Therefore, there exists a periodic C 2,Q! (f2j) solution pi of (12.31) such 
that Ui(t, x) — >■ Pi(x) as t — ?• +oo, uniformly in x G fij. 

Let now be any classical bounded solution of ( 12 .3p and let us prove that <^ > p"j in f^. 
By definition of the set is one of its connected components, and any of its connected 
components is of the type u;, + k for some k G L{L x ■ • • x L^Z. Two cases may then occur: 
either Ui is bounded, or is unbounded. 

Case 1. Consider first the case when w f is bounded. Since ^ > in (c fi$), q. t = 
on 9o;j (c dfij) and 0) = 0, it follows from Hopf lemma and the compactness of dui that 

max — — (x) < 0, 
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where v denotes the outward unit normal on dVt. On the other hand, the principal eigen- 
function ^ of ( 12. 4p is (at least) of class C 1 ^) and <pi = on dui. Hence, the quantity 

e* := sup {e G (0,e ], e^i<qi in c^"} 

is a positive real number, belonging to the interval (0, Eq\. Furthermore, E*(f>i < in UJl. 
Since e*(fi is a strict subsolution in a;, C Oj, in the sense of ( 12.51) . the strong maximum 
principle and the Hopf lemma imply that £*& < q% in Ui and 

dqt 

— < e — — on dui. 
av av 

Therefore, there exists rj > such that (e* + 77) ^ < q^ in cJJ for all 77 G [0, 770]. The definition 
of e* then yields e* = e , whence e (fi < cfi in ZDJ. The same argument can be repeated in Ui + k 
for all k G L{L x • • • x L^Z. Therefore, < qfi in f^. By comparing ^ with the solution U{ 
of the Cauchy problem (12. 6p . it follows then as in the proof of Proposition 11.11 that 

Pi < q { in fij. (2.7) 

Case 2. Consider now the case when is unbounded. For all y G C0i and R > 0, define 



'■y 



R = {z G Wj, dn(y,z) < R}, 



where <in denotes the geodesic distance inside Q, and set 



AV4> ■ V0 - C0 2 

Ai, Wiyii ,D= min — ^ . (2.8) 

*eH&K Sl .R)\{o} 



Actually, Ai M H) £> is the smallest eigenvalue of the operator —V • (AV) — C m with 
zero Dirichlet boundary condition (that is, in the Hq^^r) sense), but, since du^R may 
not be smooth in general, the eigenvalue Ai jaJi R) d may not be associated with C 1 (cJ7^r) 
eigenfunctions. We first claim that 



limsup f sup Ai MiS]J j,r>) < 0. 



To do so, let p : R — > [0, 1] be a C°°(R) function such that p = 1 on (—00, —1] and p = 
on [0, +00) and, for all y G 0Ji and i? > 0, denote 

Pv,r( x ) = P(^( x ' ?/) ~ R ) for a11 x G 
These functions are then is W l, °°(u)i). For every y G and i? > 1, the restriction of 
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the function (fiP y , R to u) it y )R belongs to Hq (u iy)R )\{0}, whence 



Ai,u, ii9ii j,£> < 



2 



/ pyfiAVlfi ■ V(tpiPy >R ) - C(<PiPv,R. 

J UJi.,. f? 



< <"^i,y,R |_ 



2 



where 



M = (l + ||V&|U~ (n .)) x _ max (A(x)£ • 

xeo, |^|=i,|^|=i 



is a positive constant which does not depend on y or R, and denotes the 

Lebesgue measure of uji^^uji^^-i. By integrating by parts, it follows then from (12. 4 p that 

Ai, WiiJ/ijKj £) S Ai,n <r D + - 



^i,y,R 



Since ^ is periodic and positive in Q{ (and then uniformly away from in each non-empty 
set of the type 

ojf : = {x G Ui, d(x,dui) > <5} 
with 5 > 0) and since Q (and hence cu;) has a smooth boundary, it follows that 



liminf ( inf \uj i>y ^i\ 1 / {<fip y , R ) 2 ) > 



while limsup K ^. +0O (sup^. |o; i)2/ijR _i| V^-rV^k-iI) = 0. Remember that A iA ,D < 0. 
Therefore, there exists Ro > 1 such that 

V/? > i? , G Wi , < (2.9) 

Let now 5 > be any positive constant such that ojf ^ and let us show that inf^g; > 0. 
Assume not and let Si > be such that F(x, s) > ^(a;)s+(Ai ( n i) £»/2)s for all (x, s) G fix [0, ei\. 
There is then a sequence (x n ) n£ N in such that qi(x n ) — > as n — > +oo. Since ^ > in cu^ 
and gj = on dwj, it follows from Harnack inequality that 

max cjj < Si for some no G N large enough. 

In particular, 

- V • (A(x)Vq t ) - C(z)ft > ^£ 5. i n ^ no ,« . (2.10) 
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On the other hand, from (12. 8p and (12. 9p . and owing to the definition of if^c^ ^ ^), there 
is G C\ (^i,x n() ,B € )\{0} (with a compact support which is included in u i>Xn ! r ) such that 



/ 



AV(f) ■ V0 - C0 2 



2 2 

Now, let a;' be any bounded open set of class C 2,a , containing the support of 0, and such 
that u' C ^i,x no ,Ro ■ ^ follows that Xi^^d — R[4>] < ^i,n 4 ,D/2. There is then a nonnegative 
and nontrivial function ip' G C 2 ' a {u') solving 

- V ■ (A(x)W) - C(*V = A w ,zV < in u' (2.11) 

with p' = on du'. Notice that <p' may not be positive in u' since u' may not be connected. 
But ip' is positive at least in one connected component u" of u' . Since min^ft > and 
ip' = on du", it follows from (I2.10p . (12.1 ip and the strong maximum principle that sip' < ft 
in a;" for all e > 0, which is clearly impossible. One has then reached a contradiction. Hence 
there holds 

inf ft > for all 5 > such that wj ^ 0. (2.12) 

a/ 

It follows then from (I2.12p . together with the Hopf lemma and the global smoothness 
of dui, that sup aw . ^ < 0. Therefore, the quantity 

e* := sup {e G (0,e ], epi < ft in u~} 

is a positive real number. From (I2.5P and the strong maximum principle, there holds e*pi < ft 
in Ui. Furthermore, we claim that 

inf (ft - e*pi) > for all 5 > such that uf ^ 0. (2.13) 

Assume not. Then there exist 8 > such that uf ^ and a sequence (y n ) ne N in uf such 
that qi(y n ) - e*pi(y n ) -)> as n -)> +oo. Write y n = + where y' n G LiZ x • ■ ■ x L^Z 
and y" G C. Notice in particular that d(y", dQ) = d(y", dfli) > 5. Up to extraction of a 
subsequence, one can assume that y" — > y^ G f2j as n — >■ +oo with 

d(yoo, <9fi) = d{y 00 , dtti) > 5, 

and that the functions x i->- ft (a; + y^) defined in f2j converge in C^ c (f2j) to a solution ft of 

-V • (A(x)Vft) = F(x,ft(x)) in H~ 

such that ft > e*pt in B(y 00 ,S) C fij with equality at y^. The strong maximum principle 
and (12.51) lead to a contradiction. Thus, the claim (I2.13P holds. As above, it follows then from 
Hopf lemma and the global smoothness of du^ that sup 3^. 9 ^ i ~ d l — < and that there exists 
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r] > such that (e* + rj)<pi < in UI for all rj G [0, rj }. Therefore, e* = e , whence e (fi < qi 
in uTi and then in Qi by repeating the argument in u>i + k for all k in L{L x • • • x L^T,. Finally, 
by comparing q,i with the solution Ui of the Cauchy problem (12.61) . the conclusion (12. 7p follows. 
Conclusion of the proof. Define the function p in Q by 



The function p is periodic, of class C 2,a (Q), and it solves (I1.13p . Furthermore, it follows from 
the previous steps that any bounded solution q of (jl,13p is such that q > p in Q. Lastly, 
let Uq : Q — > [0, M] be any uniformly continuous function such that u ^ in Q, let u be the 
solution of the Cauchy problem (11 . 14j) and let wbea connected component of Q intersecting 
the support of u . We shall prove that 



for any compact set K C uJ. Since < u(t, •) (< M) in Q for all t > and p = in Qi 
for all i ^ ij, it is sufficient to consider the case when u = + k for some i 6 L and 
some k G LxZ x • • • x L^Z. 

If is bounded, then u(l, •) > in u and maxg w < from the strong parabolic 

maximum principle. Therefore, u(l, •) > Sif>i in to for some e G (0,£o] and 



where v is the solution of the Cauchy problem ( 12.61) in ZJ with initial datum eipi in u and 
zero Dirichlet boundary condition on du. Owing to (12. 5p . v(t,x) is increasing with respect 
to t (and bounded from above by the constant M), and it converges as t — > +oo uniformly 
in uJ to a solution w of ( 12. 3p in oJ such that u> > in ZJ (whence w > in ui) and w = 
on 9a;. It follows as in the study of case 1 above that w > p in uJ, which yields ( 12. 14[) . 

Consider now the case when u is unbounded. Without loss of generality, up to a trans- 
lation of the origin, one can assume that k = and u = Ui. Choose any point yo in u 
and, from (12. 9p . let Ro > be such that Xx,u iyo Ro ,d < 0. As above, there is then a function 
(j) G Cl(u ijy0tRo )\{0} such that 



and, if u' is any bounded open set of class C 2,a containing the support of and such 
that ui' C uJi t y 0t n C co, there holds \i tC j>,D < R'[4>\ < 0. There is then a nonnegative and 
nontrivial function ip' G C 2,a (o; / ) such that 




Pi in all the sets Qi with jeL, 
in all the sets f2j with i £ I_. 




(2.14) 



w(t + 1, x) > f (t, x) for all (t, x) G (0, +oo) x w, 




V • (A(x)V<p') - ((^V 9 ' = ^i,u',d<p' in a;' 



17 



with ip' = on du' . Therefore, the function e'ip' is a subsolution of (12. 3p in u' for e' > 
small enough and one can also assume without loss of generality that e'ip' < u(l,-) in the 
compact set u' C w. Thus, there holds u(t + l,x) > v(t,x) for all (t,x) G (0, +oo) x uJ, 
where v is the solution of the Cauchy problem (12.61) in u with initial datum vq = e'tp' in u' 
and v o = in UJ\oj', and zero Dirichlet boundary condition on du. But v(t,x) is increasing 
with respect to t and bounded from above by M. It converges locally uniformly in u to a 
solution w of (12.31) in uJ, such that u> > v o in cU (whence w > in w from the strong maximum 
principle). One concludes as in case 2 above that w > p in u, which leads to ( 12. 14ft . 

Lastly, observe that, if Uq < p in Q, then u(t, •) < p in Q for all t > 0. Hence, (I2.14p 
implies that w(t, x) — >■ as t — > +oo uniformly in any compact subset K Gu, where u is 
any connected component of Q intersecting the support of uq. The proof of Theorem 11.21 is 
thereby complete. □ 



3 Relationship between the problems (11.71) with non- 
linearities f n and the problem (11.181) 

This section is devoted to the proof of Theorem 11.31 By using variational arguments, H 1 
a priori estimates and Rellich's theorem, we prove the monotonicity and the convergence 
of the principal periodic eigenvalues of the linearized operators in M. N associated with the 
functions /„, to that of problem (11. lip with zero Dirichlet boundary condition on dQ. Then, 
we show the monotonicity and the convergence of the functions p n to a solution p^ > p 
of (11.181) . The minimality of each solution p n and of p will also be used. 

Proof of Theorem 11.31 Let, for each n 6 N, Ai >n and (p n be the principal eigenvalue and 
periodic eigenfunction solving fll . lTj) . Let X^d and cp solve (ll.lip . where one can always 
assume that p > in each fi, with Xi ; n u D = ^i,d, that is i G l m %n- Call 

H^ r {R N ) = {<P G Hl 0C {R N ), <P is periodic}, L^ er (R JV ) = {</> G Lf 0C (R N ), <\> is periodic} 

and C = [0, Li] x • • • x [0, L N ]. For each nGN, there holds 

Xi t n= min R n {(f)] = R n {<p n ], 

<PeH^ N )\{0} 



where 



Rn\ 



AV(f) ■ V0 - Cn0 2 
Co 



9 

'C 

Since the sequence (Cn(x))neN is nonincreasing for each x G M N , it follows that the se- 
quence (Ai jn ) ng N is nondecreasing. 

We now claim that Ai i7l < Xi t p for each n G N. The proof is based on some standard 
comparison arguments, used in [7J |9]. We just sketch it here for the sake of completeness. 



18 



Assume that Ai >n > Ai^ for some n G N. Pick any index % G I m in- Since ( = ( n in C f2, 
there holds 

-V • (A(x)V</? n ) - COVn = Ai 5 „v9 n > Ai >D <^n in fij 

and min^-y9 n > 0. In other words, the periodic function p n is a supersolution of the linear 
equation satisfied by the periodic function ip in Since p = on <9f2j and p is (at least) 
of class C 2 (f2j), it follows from the strong elliptic maximum principle that the quantity 

e* = sup {e G (0, +oo), ep < p n in 

is actually equal to +oo. This is a contradiction since ip is positive in Therefore, 
Ai, n < Ai t £) for all n G N. 

As a consequence, the sequence (Ai jn ) ne N converges monotonically to a real number Ai >00 
such that Ai.oo < Au> Let us now show that Ai j00 = Ai^. Normalize here the eigenfunc- 
tions (p n so that ||<£>n|U 2 (Co) = 1- ^ follows that 

AVp n ■ Vp n = A i ri + / ( n pl < A li00 + / ( pl < Ai j00 + max Co- 
Co Jco Jco 

Thus, the sequence (p n ) n eN is bounded in H 1 {Cq). There exists then a function p^ 
in Hp er (R N ) such that, up to extraction of a subsequence, p> n — > weakly in ^/"^(IR^) and 
strongly in L 2 (M.^). In particular, yjoo > a.e. in M. N and ||<^oo||l 2 (c ) = 1- Let ^ be any 
compact set such that K C (IR^Q) n Co- For all neN, one has 

-(maxCn) / V>1<- (n¥>l = Ai, n - / AV(p n ■ S7<p n + / (n^l 
V K ' J K J K J Co JCo\K 

< Ai j00 + / ( p 2 n < Ai )tX3 + sup | Co |, 

Jc n \K Co\K 



whence ||v3n||L 2 (ii:) - > as n — > +oo from f ll . 16j) . Thus, p^ = a.e. in K, and then a.e. 
in M^fi and ||y3oo||L 2 (nnc ) = 1- Furthermore, since p^ G Hp er (R N ), one gets that the 
restriction of p^ to Q belongs to Hq (ft), that is the space of periodic Hl oc (Q) functions 
whose trace is equal to on dQ. Lastly, observe that 



AVp n ■ Vp n < / AVp n • Vp n = X 1>n + / ( n pl < X 1>00 + / Co<^n, 

QnCo J Co J Co J Co 



while 



Co^ 2 ,, ->■ / Co<^L = / Cv4 

Co -'Co Jcinc 



as n — > +oo. Therefore, 



AW poo ■ Vy?oo < liminf / AVp n ■ Vp n < Ai j00 + / C^ 2 



nnc n->+oo jQnCo Jnnc 



oo' 
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that is i?oo[v?oo] < Ai )(X> < X^u, where the functional i?oo is defined by 



nnc 



<P 2 
cinCo 



for all G ifg iper (f2)\{0}. But since mm^ eH i (n)\{o} Rooty] — ^i,d, one concludes 
that Ai )00 = Xi t D- In other words, Ai jTl -> Ai,d as n — >• +00. 

In the sequel, assume now that Ai,.d < 0. Consequently, for each n G N, one has 
Ai, n < Ai^ < and, from Proposition 11.11 there exists a minimal periodic solution p, n of 

-V • (A(x)Vj? n ) = /„(a;,p„) in 

< p n < M in R^. 

Fix any two integers n < m. Since 

-V • (A(x)Vp n ) - f m {x,p n ) = f n {x,p n ) - f m {x,p n ) > ill R N , 

the function p n is a supersolution for the equation satisfied by p m . From the proof of 
Proposition 11.11 there exists e m > such that all functions eip m with e G (0, e m ] are subso- 
lutions of (11.71) with the nonlinearity f m . Since mm R Np n > 0, there exists e G (0, e m ] such 
that E(p m < p n in R^. Hence, the maximum principle implies that 



v(t,x) < p n {x) for all (t,x) Gl + x 



where v is the solution of the Cauchy problem (11.51) with the nonlinearity f m and initial 
datum v = e<f m . But v is nondecreasing in t and converges as t — > +oo to a solution q 
of (11.71) with nonlinearity f m , such that < q < p n in R^. By minimality of p m (from 
Proposition ll.il) . one gets that p m < q in M. N , whence 

Pm < Pn in R N . 

In other words, the sequence of functions (p n )neN is nonincreasing and then converges point- 
wise to a periodic function Poo(x) ranging in [0, M]. 

Let us now show that poo = in M Ar \fi. By multiplying by p n the equation fll.7p with 
the nonlinearity f n , that is —V ■ (A(x)Vp n ) = f n (x,p n ), and by integrating by parts over 
the cell Co, it follows that 

/ AVpn ■ Vp n = / fn(x,p n )Pn< / f (x, p n ) p n < M X HiaX |/ |, 
J C J Co J Co U N x[0,M] 

whence the sequence (p n )neN is bounded in Hp er (M. ). Since it converges monotonically 
to Poo, one infers that p^ G -ff^^R^) and p n —> Poo as n — > +oo weakly in Hp er {R N ) and 
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strongly in L^ er (M. N ). For any compact set K such that K C (M 7V \f2) fl C , one has 



-( max g n ) I p n < - g n {.x } p n )p n = - f n {.x } p n )p n 

\Kx[0,M] JJ K J K J K 

= ~ AVp n -Vp n + / fn(x } p n )Pn 
JCo Jc \K 

< / fo(x,p n )Pn < M x max |/ |. 

Jco\K R N x[0,M] 

The assumption (j 1.1 6ft yields maxx x [ 0i M] 9n ~ > ~°° as n ~ ► +00, whence p^, = a.e. in any 
such compact K. Finally, p^ = a.e. in R-^fi. Therefore, the restriction of p^ on fl is 
in Hq (&,). Furthermore, since 

- V • (A(x)Vp n ) = F(x, p n ) in H, (3.1) 

the function is a solution of the same equation in fl in the weak HQ per (fl) sense. The 
elliptic regularity theory then implies that, up a negligible set, p^ is actually a C 2 ' a (Q) 
solution of (I1.18P and the convergence p n — > p^ holds at least in the C? oc (Q) sense. 

Lastly, let us show that > p in Q. Since p^ is nonnegative and p = in all Qi 
with i ^ one only needs to prove that p n > p in fli for all i G /_ and for all nel For 
any neN and « G L, observe that the function p n is a supersolution of C I2 . 31) in because 
it solves (13. ip in Qi and p n > on Since minQ-p n > min R ivp n > 0, there is e > such 
that e(fi < p n in fli, where <pi solves (12. 4p . Since Xi^d < 0, one can even assume without 
loss of generality that e^pi is a subsolution of (I2.3p . in the sense of (12.51) . Therefore, 

w(t, x) < p n (x) for all (t, x) 6 R + x fij, 

where iy denotes the solution of the Cauchy problem (12.61) in fli with initial datum e^. 
Since if is nondecreasing in t, it converges as t — > +00 to a solution Woo of (12.31) such that 
< e(f>i < Woo < p n in fli. From the construction of p in Theorem 11.21 and its minimality, 
one infers that p < in fij, whence 

p < p n in IV 

As a conclusion, p < p n in for all n G N, whence p < in f2. The proof of Theorem 11.31 
is thereby complete. □ 



Remark 3.1 We first show in this remark that, if F fulfills the KPP condition (II. 4p in f2, 
that is if 

F(x, s) 

s 1 — y — is decreasing in s > for all i6(l, (3.2) 



then poo = p in fl. Consider first i G /_ and let us prove that the function p solving (12. 3p 
in fli is unique. The proof is similar to the ones used for instance in [2], [5J LZ] and it is just 
sketched. Let q be any periodic solution of (I2.3P in fli. From the proof of Theorem II .2\ one 
knows that q > p in fli. But eq < p in fli for e > small enough, from the Hopf lemma 
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applied to p. Therefore, the quantity e* = sup {e > 0, eq < p in fij} is a positive real 
number. If e* < 1, then 

-V • {A(x)V{e*q}) - F(x,e*q) < in Sl u 

from (13.21) . The strong maximum principle and Hopf lemma then imply that e*q < p in fij 
and even (e* + r/)g < p in Oj for all rj G [0, r/o] and for some 770 > 0. This contradicts the 
maximality of e*. Consequently, e* > 1, whence q < p in fij and finally q = p in Actually, 
with the same arguments as those used in the proof of Theorem ll.2[ the same conclusion 
holds even if q is not assumed to be periodic. Now, if i J_, then we prove that there does 
not exist any solution q of (12 .3p that is positive in f2j (or in any of its connected components). 
Indeed, since F(x, s) < ((x)s for all and s > from (I3.2p . there holds 

-V ■ (A(x)V(elpi)) - F(x,e<pi) > -eV ■ (A(x)V&) - C(^)^; = Ai,^^ > in Q { 

for all e > 0, where ^ solves (12. 4p in O i; with Xi^d > 0. In other words, etpi is a strict 
supersolution of (12. 3p for all £ > 0. It follows with the same arguments as above or as in 
the proof of Theorem 11.21 that q < efti for all e > 0, for any solution q of (12. 3p . Therefore, 
a positive periodic solution of (12.31) cannot exist, which implies that poo = p = in Qi for 
all i I-. As a conclusion, the condition (13. 2 p implies that 

Poo = p in Vt. (3.3) 

On the other hand, we can construct examples for which (I3.3P does not hold. It is indeed 
possible to construct a situation for which A^d < and there exist an index j G {1, . . . , m} 
and so € (0,M) such that F(x,s) = Xs + s 2 for all x G fl,- and s G [0, s Q ], where A > 
denotes the principal periodic eigenvalue of the operator —V ■ (v4(x)V) in Qj with zero 
Dirichlet boundary condition on dQj. Thus, X^q^d = and j J_. Let ^ be the principal 
periodic eigenfunction of (12.41) in IX, with ( = X in Qj, such that max-^-ifj = 1. For any 
e G (0, s ], there holds 

-V ■ (A{x)V(e(pj)) - F(x,e^j) = -eV • (A(a;)V^) - Ae^- - e 1 ^) = -e 2 fi* < in 

As above, it follows from the strong maximum principle that eipj < p n in Qj for all e G (0, So] 
and for all n G N. In particular, < So^- < Poo m Qj, whereas p = in IX,- by definition. 

4 Pulsating travelling fronts and limiting minimal 

speed 

In this section, we give the proof of Theorem 11.41 We establish the relationship between the 
pulsating travelling fronts for the problems (II .ip in ~R N and (II. 8p in Q when the nonlinearity F 
is approximated with nonlinearities /„ which are very negative in IR^Q, in the sense of fl!.16j) . 
We also prove that the minimal speeds of the fronts in M> N converge monotonically to a 
quantity which is equal to in a direction e when the connected components of Q are 
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bounded with respect to e. We use especially some bounds for the minimal speeds, which 
involve some linear eigenvalue problems. 

Throughout this section, we assume that Ai^ < and e is any given unit vector of M. N . 
The functions F and f n are assumed to fulfill (I1.9P and ( I1.16p . For each n G N, one 
has Ai >n < from Theorem 11.31 The functions p n denote the minimal solutions of f ll.7j) with 
the nonlinearities f n , given by Proposition II .![ and the speeds c*(e) > denote the minimal 
speeds of pulsating fronts 4> n (x ■ e — ct, x) connecting to p n for problems (II. ip in M. N with 
the nonlinearities /„. 

Proof of part a) of Theorem 11.41 Fix any two integers n < m and let us show 
that c* m (e) < c*(e). First, remember that < p m < p n < M (from Proposition 11.11 and 
Theorem 11.31) and that both functions p m and p n are periodic in K. . Let r\ > be such 
that < r] < min R iv p m and u : M — )■ [0, M] be defined by 



u (x) 



if x ■ e > 0, 
7] if x ■ e < 0. 



Let v n and v m denote the solutions of the Cauchy problems (II. 5p with initial datum u and 
nonlinearities /„ and f m respectively. Since f m < f n , the maximum principle yields 

< v m (t, x) < v n (t, x) < M for all t > and x E R N . 

On the other hand, it follows from the results of Weinberger [3-jj that 

Vc < c^(e), sup \v m (t,x) — p m (x)\ — > as t — )■ +oo, 

x£R N ,x-e<ct 

while 

Vc > c*(e), sup v n (t, x) — > as t — )■ +oo. 

xeR N ,x-e>ct 

One infers that c* m (e) < c*(e). Consequently, the sequence (c*(e)) n6 N is nonincreasing and 
it converges to a real number c*(e) > 0. 

From the assumptions ( II. 9p and f l 1 . 1 6 j) and the regularity of F and /„, there exist a 
function F : (x,u) i-> F(x,u) and a sequence of functions (f n ) ne ^ such that: i) the func- 
tion F is defined and continuous in Q x M + , of class C°' a with respect to x G £1 locally 
uniformly in u G R + , of class C 1 with respect to m with ( := §^(-,0) G C°' a (f2), periodic 
with respect to x G f2 and F satisfies (11.91) ; m) each function / n is defined and continuous 
in l^xl,, of class C°' a with respect to x G M. N locally uniformly in u G M+, of class C 1 with 
respect to u with £ n := ^f(">0) € C°' a (M JV ), periodic with respect to x G and f n satis- 
fies (11.31) ; zzz) the functions f n satisfy A1.16H with F instead of F and g n (x,u) = f n (x,u)/u 
if u > 0, 0) = Cn( x ); ^) the function F satisfies 

F(x, u) < F(x, u) for all (x, u) G H x M + 

and F(x,u)/u is decreasing with respect to w > for all x G fi; v ) the functions / n satisfy 

fn{x,u) < J n {x,u) for all (x,u) eR N xR + 
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and f n (x,u)/u is decreasing with respect to u > for all x G WL N . 

Let Ai >n and Ai d be the principal periodic eigenvalues of problems (I1.17P and (11. lip with 
coefficients ( n and ( instead of ( n and (, respectively. Since ( n > ( n in ~R N and ( > ( in Q, 
there holds 

Ai, n < Ai 5 „ and A^d < Aix>, 

while Ai 5 „ < A^d and Ai jTl — >■ A^d as n — >• +00 monotonically, from Theorem 11.31 In 
particular, Ai /ra < Ai^ < for all n G N. Let p n be the minimal periodic solution of (II .7p 
with the nonlinearity f n , given by Proposition ll.il Actually, the function p n is unique from 
property v ) above and from [5], and it is such that p n > p n in M. N since f n > f n in M. N x M + , 
from the proof of Theorem 11.31 Let c*(e) > be the minimal speed of pulsating travelling 
fronts <p n {x ■ e — ct,x) connecting to p n for problem (11. ip with the nonlinearity f n , that 
is <j) n is periodic with respect to x G WL N , < <p n (s,x) < p n (x) and n (— 00,2) = p n (x), 
<f) n (+oc,x) = 0. As in the beginning of the proof of this theorem, there holds 

0<<(e)<c;(e), (4.1) 

since /„ < f n . Furthermore, it follows from [HI [35] that c*(e) is given by 

c;(e) = min^-, (4.2) 

A>0 A 

where k e ,\, n denotes the principal periodic eigenvalue of the operator 

Cexn ■= -V • (AV) + 2A Ae ■ V + A V • (Ae) - A 2 Ae ■ e - ( n in R N . 

Let us now show that, for every A G M, one has fc e ,A,n —> k e xo as n — > +00, where fc e ,A,D 
is the principal periodic eigenvalue of the operator 

£ e ,\,n ■= -V • (AV) + 2A Ae ■ V + A V • (Ae) - A 2 Ae ■ e - ( in Q 

with zero Dirichlet boundary condition on dQ. The proof starts as in the proof of the 
convergence Ai n — > Ai ^ in Theorem 11.31 First, it follows as in the proof of Theorem 11.31 
that the sequence (/c e ,A,ri)neN is nondecreasing and that k e ^\ tn < k e ^D for all n G M. Let Tp n 
be a principal periodic eigenfunction of £ ej A,n, that is 

Up to normalization, one can assume that H^JIi^cb) = 1- By multiplying the above equation 
by 7p n , by integrating by parts over C and by using Young's inequality, it follows that the 
sequence (jp n ) n £N is bounded in Hp er (R N ). Up to extraction of a subsequence, it converges 
weakly in Hp er (M N ) and strongly in Lp er (M. N ) to a nonnegative function Tp^ G H^ er (M. N ) such 

that H^oo ||z,2(Oo) = 1- Furthermore, since the sequence (/c e ,A,n)neN is bounded and ( n —> — oo 
as n — > +00 locally uniformly in R one infers as in the proof of Theorem 1 1 . 31 that Tp^ = 
a.e. in M 7V \fi. The restriction of Tp^ to Q is then a C 2,a (Q) periodic function such that 

£ e ,A,f#oo = ^e,A,oo^oo in Q with Tp^ = on dtt, 
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where lirrin^+oo k eXn = k eXoo < k e \ D . Since the function Tp^ is periodic, nonnegative and 
nontrivial, it follows that it is positive in f2j for some i G {1, . . . , m}, that is & e> A.oo is equal 
to the principal periodic eigenvalue fc e ,AA,-D of the operator C e ^vi in with zero Dirichlet 
boundary condition on <9f2j. But since fc e ,A,fij,-D 

> k e x,d (> ^ e ,A,oo), one concludes eventually 

that k eXoo = k e ^x,D, that is 

k e ,x,n -> k eXD asn-^ +00. (4.3) 

Assume now that all connected components of Q are bounded in the direction e, in the 
sense of ffl~T9l) . Let us show that c*(e) = 0. First, it follows from (ED, fl4~2D and (O]) that 

< c*(e) < inf ~ ke ' x ' D . (4.4) 

_ A>0 A 

On the other hand, for every A > 0, there is an index i G {1, . . . , m}, which may depend 
on A, such that k eX D = k eX Q. D and thus there is a periodic function defined in such 
that C eX o.<p = k eX DLp in fij with ip > in fi« and = on <9f2j. The function ?/> = e _A ^' e - 1 ^ 
satisfies 

- V • (A(x) W) - CfaW = m ^ (4.5) 

with > in fij and ^ = on <9f2j. Let C be any connected component of Qj, that 
is C = + for some k G LxZ x • ■ ■ x L^Z. The function t/> is positive and bounded 
in C because of (II. 19ft and since ip is bounded. It follows then from Hopf lemma and the 
smoothness of DC that there exist r > and a sequence (x n ) nS N in C such that B(x n , r) C C 
for all n G N and ip(x n ) — >■ sup c -?/> as n — >■ +00. By using the standard elliptic estimates 
and passing to the limit in (14. 5 j) in B(x n ,r), up to extraction of a subsequence, one infers 
that fc e ,A,D > liminf n ^ +00 — C(x n ) > — max rrC Finally, fc e ,A,D > — max nC f° r all A > 0, 
whence c*(e) = from (14. 4p . 

Proof of part b) of Theorem 11.41 Let c be any positive real number such that c > c*(e) 
and let (c n ) nS N be any sequence such that c n — > c as n — > +00 and c n > c* (e) for all n G N. 
Let 

u n (t,x) = 4> n (x ■ e - c n t,x) 
be pulsating travelling fronts for (II. ip in R^ with nonlinearity f n , such that 



= n (+oo, x) < 0„(s, x) < 4> n (— oo, x) = p„(x) < M for all (s, x) G R x 



Actually, from [TB], each solution w n satisfies (n„)( > in R x M N . 

On the one hand, since < u n (t, x) < p n (x) in R x M. N , Theorem 11.31 implies that u n — > 
in L^ 0C (R x (R Af \fi)). On the other hand, since f n (x, s) = F(x, s) for all (x, s) G H x R + , it 
follows from standard parabolic estimates that there exists a function u : R x Q — y [0, M] 
such that, up to extraction of a subsequence, u n — > u as n — > +00 in C£ and C% in R x £7, 
where u obeys 

tt t - V ■ (A(x) Vu) = F(x, tt) in R x Q 

and < u(t,x) < Poo{x) < M for all (t, x) Glxfl, under the notation of Theorem 11.31 In 
particular, the function u can be extended continuously by on R x dVt and, from parabolic 
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regularity, the function u is a classical solution of (11. 8p in R x f2 (of course, one could also 
extend u by in R x (R Ar \fi) and u would then be continuous in R x M N ). Moreover, the 
equalities 



u n (t H , x) = u n (t, x — k) in R x 

V c n J 



DiV 



for all k G L{L x • • • x carry over at the limit, whence u(t + (k ■ e)/c, x) = u(t, x — k) 
in R x Q for all k G LiZ x • • ■ x L^Z. In other words, the function u can be written 
as u(t, x) = (fi(x ■ e — ct,x) in R x Q where (ft : R x Q — > [0, M] is such that (p(s, •) is 
periodic in Q for all s G R. Lastly, since all functions u n are increasing in time in R x M. N , 
the function u is such that n ( > in R x fi. From the previous observations and parabolic 
regularity, there are then two periodic functions u ± defined in Vt such that < u~ < u + < 
in Q, u(t,x) — > m ± (x) as t — > ±oo in Cf oc {VL) and u^ 1 obey 

-V • {A{x)Vu ± ) = F(x, tt ± ) inH, 
u ± = on <9f2. 

Let now any index % G that is Ai^^d < in the sense of (11.121) . From the proof 
of Theorem 11.21 there is a minimal periodic solution of (12. 3p . Furthermore, in there 
holds pi = p < Poo < p n for all n G N, under the notation of Theorem 11.31 Therefore, one 
can always shift in time the functions u n so that, say, 

u n (0, x) dx = - / Pi{x) dx, 

where we recall that C = [0, L x ] x ■ • • x [0, L N ]. From Lebesgue's dominated convergence 
theorem, the function u satisfies the same equality at the limit, whence 

0< / u~(x) dx < - / Pi{x) dx < / u + (x)dx 



Jc nQi 2 jc nQi Jc nQi 

by monotonicity of u with respect to t. The minimality of pi and the strong maximum prin- 
ciple imply that u~ = in while u + > in Q iy again from the strong maximum principle. 
If we further assume that F satisfies the KPP assumption (13.21) in Q (or just in Qi), then it 
follows from Remark 13. II that the solution of ( 12. 3p is actually unique, whence u + = Pi in Qi 
in this case. 

Proof of part c) of Theorem 11.41 Firstly, it follows from [31 [35] that, for each n G N, 

— k x — k x 

for some A n > 0, where k e ^ n denotes the principal periodic eigenvalue of the operator 

Cexn ■= -V • (AV) + 2A Ae ■ V + A V • (Ae) - A 2 Ae ■ e - ( n in R N . (4.6) 

Since, as above, fc e ,A,n — > k e ^D as n — > +oo nondecreasingly for every A G R, where /c e ,A,D is 
the principal periodic eigenvalue of the operator 

£ e ,\,n ■= -V • (AV) + 2A Ae ■ V + A V ■ (Ae) - A 2 Ae ■ e - C in 
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with zero Dirichlet boundary condition on <9f2, it follows that 

X n \ n a>o A 

for all n £ N, whence 

c *( e )> inf Z^. (4.7) 

A>0 A 

Furthermore, the maps A h-> — & e ,A,n are all convex and their derivatives at A = are all equal 
to 0, see [21 E]. In particular, for every nGN, — & e ,A,n is nondecreasing with respect to A > 
and — k e ,x, n > —k e ,o,n = — Ai,„ for all A £ M. By passing to the limit as n — > +oo pointwise 
in A, one gets that the map A i— > — k £) \,D is convex in R, nondecreasing in M + , and there 
holds —k e Xl D > —k e ,o,D = — Ai,d > for all A £ R. Notice here that, if assumption (I1.19p 
is made, then — k e% x,D — max fiC f° r an A, under the notation used in the proof of part a). 
Therefore, the infimum in ( 14. 7p is not reached in general. 

Because of (14. 7p . the inequality —k e x,D > — Ai^ > and the limit lim n _^. +t>0 & e> A,n = & e ,A,D 
for all A, it follows that, in order to show the positivity of c*(e), it is sufficient to prove that 
there exist A > and a > such that 

- K Xn >a\ 2 for all A > A and for all n £ N. (4.8) 

Of course, from the proof of part a), this cannot be always true. However, assuming from 
now on that A is constant, we shall now show that ( 14. 8 j) holds under conditions ( 11.201) 
or (ll.2ip . Assume first that there exist a unit vector e' ^ ±e and two real numbers a < b 
such that (ll.20p is fulfilled, that is 

^ 5 e ', a ,6 := {x G 1^, a < x ■ e' < b}. 



For any A > 0, let ip\ be the function defined in S e i >a) b by 

U*) = e AVe ' } cos x (x • e' - i±*)) , 

where A' = A (Ae • e')/(Ae' ■ e'). The function ^ is bounded and of class C^CS^^b), it is 
positive in S e > >a j) and vanishes on dS e ^ a ,b- Furthermore, since ( n = ( in Q D S'e'.a.&j it is 
straightforward to check that 

/ir 2 (Ae' ■ e') \ 

for all n £ N, where a = (Ae-e)/2 — (Ae-e') 2 /(2Ae'-e') > from Cauchy-Schwarz inequality, 
since the unit vectors e and e' are not parallel. Since ( is bounded in Q, it follows that there 
exists A > such that £ e ,A,nV'A < ~ a A 2 V'a i n £e',a,& for all A > A and n £ N. This 
inequality yields (j4.8j) . as in the course of the proof of Proposition 11.11 We just sketch the 
proof here. Fix any A > A and n £ N and let <p n be a principal periodic eigenfunction of the 
operator £ e ,A,n- Namely, C eXn 4> n = k et \, n (fi n and 4> n is periodic and positive in WL N . Define 

e* = sup {e > 0, eipx < 4>n in S e >, a ,b}- 



27 



Owing to the definition of ip\ and the uniform positivity of n , the quantity e* is a posi- 
tive real number. Furthermore, s*ip\ < <p n in 5 e / jttj 6 and there is a sequence (x m ) me ^ of 
points in S e i^ b such that \imm{ m ^ +OD d(x m ,dS e ^ atb ) > 0, \im. m ^ +oa (e*ip\( (x m )) =0 

and liminf m ^ +00 £ eiAjn (£:*^A - 4>n){xm) > 0. Since there holds C e ,\,n4>n{x m ) = k e ,x,n4>n{x m ) 
and £ e ,\,nip\(x m ) < —a>\ 2 ip\(x m ) for every m G N, one concludes that fc e> A,n < —a A 2 , that 
is (14. 8p . This yields the desired inequality c*(e) > 0, as already emphasized. 

Assume now that there exist a unit vector e', a point Xq G and a positive real 
number r such that e' is an eigenvector of A with e' ■ e / and (ll.2ip holds. Let (5 > be 
such that Ae' = fie! . Since the matrix A is symmetric, there is an orthonormal family of 
eigenvectors e^, . . . , e^ r _ 1 of A in such that e< • e' = for all 1 < i < N - 1. Even if it 
means decreasing r > in (ll.2ip . one can assume without loss of generality that 

tt D G e ,r ■= {x G R N , \{x - x ) ■ e'i\ < r for all 1 < i < N - l}. 



For any A > 0, let ip\ be the function defined in C e >^ r by 



M*)= II e^-Ocos(- 



Ki<Af-l 



COS r^-^)-^) - 

V 2r 

where A- = X(Ae ■ e'^/i^Ae^ ■ e'f). The function ip\ is bounded and of class C°°(C e ' )r ), it 
is positive in C e i, r and vanishes on dC e i )T . Furthermore, since ( n = ( in Q D C e / >r , it is 
straightforward to check that 

Ki<iV-l 



in C e / , 



for all neN, where a = /3(e • e') 2 /2 > since /3 > and e • e' ^ by assumption. Thus, one 
concludes as above that there is A > such that C ej x, n ipx < — a A 2 ?/>a in C e ', r for all A > A 
and n G N. This yields (14. 8 j) and finally c*(e) > 0. The proof of Theorem 11.41 is thereby 
complete. □ 

Remark 4.1 In the case when the functions /„ fulfill the KPP assumption ( 11. 4ft . then c*(e) 
is given by an explicit variational formula. Namely, under assumption ( 1 1.4ft for the func- 
tions /„, it follows from the proof of part a) of Theorem 11.41 with the choices f n = f n 
and F = F that c*(e) < inf A>0 —k^jj/X, because of (14.41) . On the other hand, the reverse 
inequality (14. 7p always holds, from the proof of part c) of Theorem 11.41 As a conclusion, the 
assumption ( II. 4p for the functions /„ yields 

*/ \ • c ~k e \n 
c (e) = mt . 

A>0 A 
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